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Necessary  conditions  are  derived  If  very  long  plasna  filaments 
(strlatlons)  are  to  reaaln  Magnetic  flux  tube  aligned  In  distorted 
geoaagnetlc  fields  subsequent  to  nuclear  burst.  Special  E,  fields 
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20.  ABSTRACT  (Continued) 

fields  cannot  exist,  long  nuclear  burst  strlatlons  In  the  magnetosphere 
will  gradually  shear  across  magnetic  flux  tubes,  break  up  Into  several 
shorter  segments,  descend  Into  the  lower  Ionosphere,  and  be  consumed 
faster  there  by  molecular  deionization. 
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An  Important  question  concerning  the  late-time  evolution  of 
nuclear  burst  plasma  plumes  and  strlatlons  Is  how  accurately  the  plasma 
filaments  stay  aligned  with  magnetic  flux  tubes  as  the  filaments  move 
about.  Because  of  various  forces  acting  on  the  plasma,  time-varying 
currents  are  induced  in  the  plasma  that  more  or  less  distort  the  magnetic 
flux  tubes.  The  existence  of  plasma  currents  also  means  that  the  Ions  and 
electrons  do  not  move  with  exactly  the  same  velocity.  The  question  we 
wish  to  address  here  Is  the  following:  given  a  very  long  plasma  filament 
aligned  perfectly  with  a  particular  magnetic  flux  tube  at  some  Initial 
time,  how  accurately  will  the  plasma  elements  of  the  filament  be  aligned 
along  a  single  flux  tube  at  later  times? 

There  Is  a  strong  tendency  for  an  Ideal  MHO  plasma  to  move  about 
In  such  a  manner  as  to  stay  aligned  with  a  magnetic  flux  tube.  For 
example,  a  low  energy  density  plasma  moves  transverse  to  B  with  a  velocity 
nearly  equal  to  cExB/B2,  and  It  has  been  shown1.2  that  one  can  consider 
the  magnetic  flux  tubes  to  move  with  this  same  velocity  If  VxE(  happens  to 
be  zero  or  parallel-to-B  at  all  times.  Within  these  approximations,  one 
can  then  consider  the  plasma  filaments  and  magnetic  flux  tubes  to  remain 
aligned  at  all  times. 

However,  the  existence  of  plasma  currents  Implies  that  the  Ions 
and  electrons  do  not  both  accurately  follow  the  cExB/Bz  drift  velocity, 
and  that  resistive  or  anomalous  Ef  fields  may  occur  whose  VxEj  may  not  be 
parallel-to-8.  Consequently,  a  long  plasma  filament  originally  aligned 
with  a  magnetic  flux  tube  will  not  necessarily  map  onto  a  single  flux  tube 
at  later  times.  It  Is  this  time-evolving  shear  of  a  plasma  filament  onto 
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several  flux  tubes  that  we  investigate  here  because  it  could  substantially 
speed  up  the  late-time  decay  of  nuclear  burst  striations  by  allowing  a 
long  filament  to  break  up  into  several  segments  that  can  more  easily  fall 
into  the  lower  Ionosphere  and  be  consumed  by  molecular  chemistry  there. 

We  show  that  special  E(  fields  would  be  required  if  striations 
are  to  remain  magnetic  flux  tube  aligned.  We  have  not  yet  been  able  to 
demonstrate  that  such  special  E,  fields  can  appear  in  reality.  Some 
observed  auroral  E(  fields  seem  similar  to  our  special  Ej  fields.  The 
observed  westward  motion  of  the  Checkmate  and  Kingfish  beta  tubes  may  be 
consistent  with  a  temporary  appearance  of  these  special  E(  fields,  but 
further  detailed  computations  are  required  to  reach  a  definite  conclu¬ 
sion. 


♦ 

If  the  above  special  Ejj  fields  cannot  develop  in  nuclear  burst 
striations,  then  It  seems  likely  that  long  striations  will  gradually  shear 
across  adjacent  flux  tubes  In  those  regions  of  the  striation  where  J.  is 
non-zero,  with  a  transverse-to-B  shearing  velocity  equal  to  cJx/eN.  The 
effect  on  the  striated  nuclear  burst  plasma  would  be  to  decrease  the  mean 
electron  density  N  and  rms  electron  density  fluctuations  oN  at  very  high 
altitudes,  to  Increase  the  N  and  at  Ionospheric  altitudes,  and  to  ac¬ 
celerate  the  overall  plasma  decay  due  to  the  more  rapid  molecular  chemis¬ 
try  at  ionospheric  heights. 

We  show  that  In  Ionospheric  plasmas  the  electrons  move  In  direc¬ 
tion  of  Jx  relative  to  the  background  neutrals.  This  Is  contrary  to  the 
folklore  of  plasma  theorists  that  "the  Ions  carry  the  transverse-to-B 
currents".  Thus,  predictions  of  Image  striations  In  the  E-layer  may  be  a 
spurious  result  of  theorists  If  the  above  special  E(  fields  could  appear 
In  actual  plasmas. 
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The  author  benefited  from  many  lengthy  and  stimulating  discus¬ 
sions  of  this  material  with  F.  Fajen,  R.  Stagat,  W.  White,  C.  Longmlre  and 
J.  Sperling.  Briefer,  but  useful,  discussions  were  also  held  with  0. 
Mai oof,  T.  Mazurek,  0.  Sowle,  L.  Wlttwer  and  N.  Krall.  The  author  Is  also 
Indebted  to  W.  Chesnut,  J.  Workman  and  C.  Prettie  for  Information  concern¬ 
ing  Checkmate  beta-tubes  and  Barium  cloud  strlatlon  behavior. 

Due  to  the  vociferous  discussions  unleashed  by  a  preliminary 
version  of  this  report  (dated  30  June  1983),  the  present  version  has  been 
revised  to  clarify  the  need  for  special  Et  fields  If  “Image  strlatlons" 
are  to  be  suppressed. 
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SECTION  1 

DRIFT  VELOCITY  OF  ELECTRONS,  IONS,  AND  FLUX  TUBES 


We  assume  the  electron  fluid  and  ion  fluid  can  be  described  by 

♦ 

the  usual  fluid  MHO  equations  (at  least  in  the  transverse-to-B  directions) 

and  Maxwell's  equations.  The  two  momentum  equations  for  the  electron 

♦ 

fluid  and  ion  fluid  can  be  inverted  to  express  the  transverse-to-B  elec- 
♦  ♦  ♦ 
tron  velocity  V  .  and  ion  velocity  V.,  in  terms  of  E,  ,  the  neutral  wind 

velocity  V^,  the  Pedersen  conductivity  a and  Hall  conductivity  ch  ,  ion 
fluid  acceleration,  gravity,  pressure  gradients,  etc.  The  exact  expres¬ 
sions  are  algebraically  very  complex,  and  are  given  on  pages  617-621  of 
Reference  3.  For  convenience  of  the  reader,  this  section  of  Reference  3 
Is  reproduced  as  an  addendun  to  this  note.  For  particular  model  problems, 
one  usually  will  simplify  the  full  equations  (12-133,  12-134,  and  12-135) 
for  the  relevant  physical  problem  being  Investigated. 

Newcomb1  has  shown  how  one  may  ascribe  a  velocity  field  V.  for 

the  motion  of  magnetic  flux  tubes  when  Ex  and  E(  fields  are  present. 

Because  the  flux  tubes  are  not  true  entities,  there  is  some  degree  of 

arbitrariness  to  this  velocity  field  which  one  can  exploit  in  particular 

♦  ♦  ♦ 

problems,  e.g.,  by  matching  to  or  on  some  convenient  surface 
that  cuts  through  all  the  magnetic  field  lines  (see  Figure  1). 

♦  ♦ 

Given  Ej_  and  Ef  over  the  magnetic  volume  of  interest,  Newcomb 
showed  that  the  flux  tube  velocity  is: 


Figure  1.  Oistorted  geomagnetic  field  lines  and  an  example  of  a  chosen 
surface  that  cuts  all  relevant  field  lines.  A  particular  field 
line  is  identified  by  the  coordinates  x§  and  y$  in  the  surface 
where  the  field  line  cuts  through  the  surface,  and  is  the 

distance  along  the  field  line  from  that  surface  point. 


Here  x,  and  y,  are  surface  coordinates  in  a  chosen  surface  that  cuts 
s  s 

through  all  the  magnetic  field  lines,  and  the  line  integral  /E,dt,  is  to 
be  evaluated  along  each  magnetic  field  line,  starting  at  the  point  where 

the  field  line  cuts  through  the  surface  (i.e.,  Jt,= 0  there).  The  arbitrar- 

♦  ■ 

iness  of  the  is  given  by  the  term  V^f^.y^A^O),  where  f(xs,y$  )  may 
be  any  desired  function  defined  on  the  chosen  surface.  Usually  it  is 
desirable  to  place  the  surface  where  E,  is  negligible  and  choose  f=0;  near 
this  surface  Equation  1  then  yields  V^cE^xB/B2,  but  more  distant  points 
from  this  surface  will  also  have  a  contribution  from  the  second  term,  of 
Equation  1. 

Equation  1  satisfies  Faraday's  law  because: 


38  ♦ 

—  =  -  cvxE 


*  -  cvx[  E  +  ^(xs,ys,A|)] 

■  Vx(VfxB)  (2) 

Newcomb  points  out  that  the  last  step  is  valid  if  we  choose 

V(WV  =  ‘  ^xs’ys’V 

which  on  integration  along  the  field  lines  starting  from  our  surface  (at 
Jt,=0)  yields 


f(xs,ys.*,)  ■  f(xs,ys,£1*0)  -  /  E^t, 
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where  f(xs,ys>£(*0)  may  be  chosen  arbitrarily  on  the  surface  that  cuts  all 
field  lines.  Consequently  from  Equations  2  and  3: 


*  -  f*  m  L  ■  \  ['  E.d*. ♦  v(w*i  ■  °>)  <5> 

Thus  this  velocity  field  describes  exactly  the  evolution  of  the  magnetic 
field  if  we  are  given  and  Ej  at  all  times  and  at  all  locations. 

How  much  a  plasma  filament  becomes  misaligned  from  a  flux  tube 
depends  on  how  much  the  ion  velocity  V..  or  electron  velocity  Val  deviates 
from  the  flux  tube  velocity  vfx  at  different  points  along  the  flux  tube, 
given  that  the  plasma  velocity  and  flux  tube  velocity  are  set  equal  at  a 
chosen  surface  that  cuts  across  all  field  lines.  Comparing  Equations 
12-133  and  12-134  of  the  addendum  with  Equation  5  above,  it  seems  likely 
that  some  deviations  exist  between  these  three  transverse-to-8  velocity 
vectors.  However,  these  general  equations  are  so  complex  that  it  is 
difficult  to  come  to  specific  conclusions.  Thus  we  will  examine  a  couple 
of  simple  problems  to  exemplify  possible  plasma  motion. 
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SECTION  2 

EXAMPLE  OF  IONOSPHERIC  PLASMA  MOTION 

Let  us  consider  a  plasma  in  the  lower  ionosphere  where  the 
dominant  forces  are  collisions  of  ions  and  electrons  with  neutrals.  Thus 
in  Equations  12-125  and  12-126  of  the  addendim  we  neglect  the  inertial 
terms,  the  pressure  gradients,  and  gravitational  force,  so  momentum 
balance  for  the  ion  and  electron  fluids  is: 

O  -  el»(M,^/c)  ♦  »,„»,(*„-?,)  -  ve)pe(V,-Ve)  (6) 

0  -  -eN(E+Ve*8/c)  ♦  Ve(Vo-y  *  ve1pe(vrve)  (7) 

Addition  of  these  two  equations  yields 

0  •  W  *  Vt<VV  +  %»e(»0-;e>  '8> 

where  we  have  employed  the  definition 

J  •  |1!  (vr»e)  (9) 

♦  ♦  ♦  ♦  ♦ 

He  can  now  express  V^,  E±  and  in  terms  of  J^,  the  neutral  wind 
velocity  V^,  and  the  magnitudes  of  the  ratios  of  collision  frequencies 
over  gyro  frequencies  (n^-v^/a^  r^-v^/a^  and  ne4=vei/ne  )  by  the 
following  procedure:  (a)  use  Equation  9  to  eliminate  Vg  in  Equation  8  and 
solve  for  V^;  (b)  use  the  resultant  equation  and  Equation  9  to  express 
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In  terms  of  and  V^;  (c)  then  use  Equations  6  or  7  to  find  E^,  and 
Equation  5  to  find  V^.  The  resultant  expressions  are: 


-  V 

+  1  . 

cJxb  .  neo 

c\ 

(10) 

04- 

n1o4T'eo 

eN  "lo^eo 

eN 

♦ 

•  V 

+  1  . 

c  Jx  b  11 1  o 

C*4. 

(11) 

04- 

^io^eo 

eN  n.„ -*n 

to  eo 

eN 

'fj-  <u- 


*-  ♦  [  — ♦  *„]  ^  -  |Mxf(xs,ys,0) 


"lo%o 


B 


n1o-neo  cJJ-  *  ♦ 


n1o+neo 


eN  B 


+  -/  Ei<“. 


(12) 


-► 

E. 


V8  A  r  1+n1oneo  .  1  BJ4-  *  n1o”T,eo  JxB 

-  l-  -  eii  eN  nl0%0  eN 


n1o+  neo 


(13) 


No  approximations  have  been  made  In  deriving  Equations  10-13  from 
Equations  5-9,  except  that  the  electron  density  Ng  Is  very  nearly 
equal  to  the  Ion  density  N^. 

♦ 

Note  that  the  first  two  terms  in  Equations  10  and  11  for  and 
V„.  are  the  same.  The  last  term  of  Equations  10  and  11  shows  the  velocity 
components  of  and  In  direction  of  because  h^o«*1000neo  at  all 
altitudes,  this  shows  that  the  electrons  carry  essentially  all  of  the 
transverse-to-8  current  density  relative  to  the  neutrals. 
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More  precisely,  the  vector  dot  products  of  Jx  with  Equations  10 

and  11  show  that  In  a  reference  frame  moving  with  the  local  neutral  wind 

velocity  VQX>  It  Is  essentially  only  the  electrons  that  move  In  direction 

of  and  the  Ions  have  essentially  no  motion  In  the  direction.  The 

writer  dearly  hopes  that  consideration  of  Equations  10  and  11  will  lay  to 

♦ 

rest  the  hoary  tale  of  plasma  theorists  that  "the  ions  carry  the  Jx  be¬ 
cause  the  ion  Pedersen  conductivity  is  a  1000-fold  larger  than  the  elec¬ 
tron  Pedersen  conductivity".  Ionospheric  theorist  seem  more  at  home  with 

equations  expressed  in  terms  of  E.  and  Exb;  such  an  alternative  derivation 

4“  *►  ♦ 

is  given  in  the  Appendix,  again  showing  that  Jx»Vix  is  very  small  compared 

to  in  the  reference  frame  of  the  neutral  fluid. 

4-  e-u 


Because  the  J,  is  also  carried  primarily  by  the  electrons,  Equa- 
*  -► 

tions  10  and  11  show  that  the  current  density  J  is  a  result  of  primarily 
electron  motion  around  the  whole  current  circuit. 

Let  us  consider  now  the  conditions  needed  such  that  a  plasma 

filament  will  stay  field-aligned  during  its  drift  motion.  This  means  that 

♦ 

we  would  like  the  flux  tube  velocity  V  f  to  be  equal  to  either  the  ion 

♦  +.-W 

velocity  or  the  electron  velocity  V^.  Comparing  the  second  term  of 

Equation  12  with  that  of  Equations  10  and  11,  we  see  that  we  must  have 

neo*°  and  nel“°*  ^*e*»  e^ectron  collision  should  be  negligible.  (This  is 

In  line  with  the  usual  expectation  of  theorists.)  In  this  limit,  the 

next-to-last  term  of  Equation  12  reduces  to  -cJx/eN,  so  the  magnetic  flux 

tubes  could  move  with  the  electrons  If  the  last  term  of  Equation  12  were 

negligible,  e.g..  If  E.  were  zero.  The  flux  tubes  would  distort  if  they 

♦  ■ 

follow  the  electron  V  ,  that  Is  consistent  with  the  plasma  motion  and  the 
♦  e-*- 

Jx.  In  this  case  (Ef -0  and  ne0*nej  *0)  "Image  strlatlons"  can  form  In 
regions  where  Is  non -zero. 


Alternatively,  the  magnetic  flux  tubes  might  follow  the  Ion 
♦  ♦ 

motion  (l.e.,  w  ).  In  addition  to  requiring  no  electron  collisions, 

the  continuing  alignment  of  the  Ion  fluid  and  a  flux  tube  requires  that 

the  last  two  terms  of  Equations  12  cancel  each  other.  Crossing  these  last 
♦ 

two  terms  with  B,  this  requires: 


Ai 

’a./'  E 


l^l 


^o^eo  #  JxB 
"lo%o  ’  *K 


(14) 


Thus  the  presence  of  current  loops  In  this  case  requires  the  appearance  of 
♦ 

special  E  fields  If  the  Ion  fluid  Is  to  stay  flux-tube  aligned.  These 

♦  I 

special  E,  fields  do  not  seem  to  be  associated  with  any  resistive  effects 

■  ♦ 

associated  with  J(  .  In  fact.  Equation  14  must  hold  in  the  limit  of  no 

electron  collisions  If  the  Ion  fluid  Is  to  stay  flux  tube  aligned.  In 

this  case  no  "Image  strlations"  can  form  because  the  electrons  move  across 

♦ 

the  magnetic  flux  tubes  In  direction  of  Jx. 


It  Is  noteworthy  that  the  r.h.s.  of  Equation  14  is  exactly  the 

same  term  as  the  last  term  of  Equation  13  for  E,.  This  last  term  for  E. 

appears  to  be  that  part  of  E.  required  so  that  the  electrons  can  "ExB" 

drift  across  the  plasma  and  thus  generate  the  current  density  Jx  ,  l.e., 

the  last  term  of  Equation  11.  Denoting  this  last  term  of  Equation  13 

♦  ♦  ♦ 

(which  Is  along  the  direction  JxB)  as  E^g  >  we  can  then  rewrite  Equation 
14  as: 


(15) 


Taking  the  derivative  of  this  equation  w.r.t.  ,  this  appears  to 
Imply: 


(16) 


7x(El  +  E0JxB) 


Thus  we  see  that  the  purpose  of  the  special  Ej  ,  required  by  Equation  14, 

is  to  cancel  any  non-zero  curl  of  the  last  term  of  E.  in  Equation  13  and 

thereby  suppress  the  associated  3B/3t.  That  is,  that  portion  of  the  total 
♦ 

E  field  that  appears  in  Equation  16  must  be  the  gradient  of  a  potential: 


Ei  +  E-U)xB 


(17) 


Although  the  above  discussion  may  strike  the  reader  as  somewhat 
bizarre,  there  appears  to  be  a  theoretical  similarity  with  Ef  fields  that 
are  observed  to  be  associated  with  auroral  arcs.  For  example,  see  the 
discussion  of  auroral  E(  fields  by  Hallinan4. 


The  time  evolution  of  the  plasma  motion  and  magnetic  field 

♦  ♦ 

configuration  can  be  computed  from  Equations  10  and  12  for  and  , 

plus  Ampere's  Law  (J-VxB/4w)  and  an  expression  for  Ej  (e.g..  Equation  14 

or  electron  momentum  balance  parallel -to-B).  The  distorting  magnetic 

♦  «► 

field  lines  are  advanced  In  time  by  the  Integral  of  VfJ_,  and  B  is  computed 

from  the  spatially  distorted  field  lines.  Note  that  Equations  11  and  13 

are  not  necessary  to  solve  for  the  plasma  and  magnetic  field  evolution, 

but  are  merely  auxiliary  equations  that  allow  the  computation  of  and 

E.  If  these  quantities  are  desired.  This  procedure  conforms  with  Fara- 
♦  ♦  ♦ 
day's  law  (3B/3t*-cVxE)  because  of  the  above  choice  for  VfJ_. 
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SECTION  3 

A  SIMPLIFIED  DETAILED  EXAMPLE 


The  rising  Kingfish  fireball  and  associated  kinked  beta-tube 
ought  to  be  described  fairly  well  by  Equations  10-13.  Thus  if  we  were 
given  the  density  pQ  and  velocity  VQ  of  the  neutrals,  and  the  ion  number 
density  N  in  the  fireball  and  the  beta- tube,  then  the  above  equations 
should  predict  the  evolution  of  the  plasma  location  and  magnetic  field 
configuration.  The  beta  tube  behavior  was  clearly  visible  below  -120  km 
altitude,  and  clearly  showed  that  field  lines  leaving  the  high  altitude 
fireball  were  highly  distorted  down  to  about  85  km  altitude  for  about 
three  minutes  after  the  burst.  A  peculiarity  of  the  data  is  that  the 
footpoints  of  these  field  lines,  at  -85  km  altitude,  moved  about  50  km 
westward  of  the  fireball  magnetic  meridian,  as  well  as  lagging  behind  the 
fireball  rise  in  their  northward  motion.  The  upward  bend  of  these  field 
lines  at  -85  km  altitude  was  as  much  as  45*,  so  a  substantial  current  loop 
was  set  up  linking  the  high  altitude  fireball  plasma  to  the  beta-tube  kink 
down  at  85  km  altitude. 

An  explicit  computation  of  the  Kingfish  behavior  is  beyond  the 
scope  of  the  present  paper,  but  we  shall  consider  a  model  problem  that  I 
believe  retains  the  essence  of  the  Kingfish  problem,  and  therefore  will 
allow  us  to  compare  with  the  observed  field  line  behavior.  Figure  2  shows 
the  model  plasma  and  neutral  wind  configuration.  The  figure  might  corres¬ 
pond  to  a  nuclear  burst  in  the  lower  ionosphere  at  the  geomagnetic  equa¬ 
tor,  wherein  the  central  plasma  represents  the  fireball  and  the  northern 
and  southern  plasmas  represent  the  northern  and  southern  conjugate  region 
(NCR  and  SCR)  beta-tube  plasmas.  Here  we  are  In  a  reference  frame  moving 
upward  with  1/2  of  the  fireball  velocity. 
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(a).  Initial  configuration.  View  from  east. 


(b).  Final  quasi-steady  state  configuration.  View  from  east. 


Figure  2.  Model  plasma  that  exemplifies  the  Klngflsh  bent  beta-tube 


In  the  model  problem  shown  In  Figure  2,  there  are  three  high 
density  plasma  elements  along  the  Initially  undistorted  magnetic  flux 
tube,  embedded  in  a  background  of  low  density  plasma.  As  an  example,  the 
three  high  density  regions  might  have  10®  ions/cm3,  and  the  surrounding 
low  density  plasma  might  have  10s  ions/cm3.  The  z-axis  Is  vertical  (e.g., 
located  at  the  geomagnetic  equator),  while  the  x-axis  is  horizontal  and 
pointing  northward.  We  assume  the  length  of  the  central  plasma  element  is 
twice  the  length  of  the  two  equally  long  northern  and  southern  elements. 
The  neutral  fluid  density  is  assumed  uniform  everywhere,  and  there  is  an 
upward  neutral  fluid  velocity  VQ  across  the  central  plasma  element, 
and  an  equal  but  downward  velocity  -V  across  the  northern  and 
southern  plasma  elements.  Consequently  there  is  an  east -west  vertical 
plane  of  symmetry  across  the  middle  of  the  central  plasma  (say  at  x=0), 
and  we  need  consider  only  the  details  of  the  northern  half  of  the  model 
problem. 

Figure  2a  shows  the  configuration  at  the  initial  instant  when 
the  vertical  neutral  winds  are  turned  on.  If  the  three  plasma  elements 
stayed  aligned  within  a  common  distorted  magnetic  flux  tube,  it  would 
eventually  take  on  the  quasi-steady  state  configuration  shown  in  Figure 
2b. 

Figure  3  shows  the  current  configuration  for  the  case  when  the 

transverse-to-8  plasma  cross-section  is  rectangular  (say  lOkmxlOkm).  The 

top  view,  shown  in  Figure  3a,  shows  a  westward  J.  across  the  central 

plasma  element;  the  associated  downward  J^xB  force  balances  the  large 

upward  ion-neutral  drag  force  (the  second  term  of  Equation  8)  and  thereby 
♦ 

determines  at  any  instant.  When  this  Jx  reaches  the  west  face  of  the 

central  plasma  element,  it  turns  into  Jj  that  runs  down  the  field  lines  to 

the  northern  (and  southern)  plasma  element.  There  the  J.  turns  into  an 

♦  ♦  ■ 

eastward  J±,  whose  associated  upward  J^xB  force  balances  the  downward 
Ion-neutral  drag  force  there.  When  this  Jx  reaches  the  east  face  of  the 


northern  plasma  element.  It  turns  Into  J(  and  runs  back  up  the  field  line 
to  the  central  plasma  element  to  complete  the  current  circuit. 

Figure  3b  shows  the  distorted  flux  tube  and  Jj_  directions  as 
seen  from  the  east.  Note  that  the  field  lines  are  bent  in  the  regions 
where  flows.  If  the  field  lines  are  bent  an  angle  $  across  the  region 
of  a  plasma  element,  one  can  relate4  the  field-line-integrated  J.  to  the 
kink  angle  0  across  the  plasma  element  via  Ampere's  Law  (J=VxB/4ir): 

/  JLdii  *  7-  sing  (18) 

1  Air 


We  can  use  this  expression  to  relate  to  the  kink  angle  0  across  each 
plasma  element  by  field-line-integrating  Equation  8  along  each  plasma 
element.  This  yields  for  each  plasma  element 


•fnioeNVodJli  "  cBsin0/4ir 
/n^eNdt, 


(19) 


Assuming  now  that  the  flux  tube  distorts  with  the  plasma  velocity  V-x,  the 
transverse-to-B  displacement  of  each  flux  tube  element  is  given  by  /V.^dt, 
and  the  resultant  kinked  flux  tube  geometry  determines  the  various  kink 
angles  &(t)  as  a  function  of  time  at  each  plasma  element.  Thus  Equation 
19  allows  a  complete  computation  of  the  flux  tube  evolution  in  time 
(assuming  that  It  moves  with  the  local  Ion  fluid  velocities  V14_).  In  par¬ 
ticular,  for  our  simple  model  problem.  Equation  19  shows  that  a  steady- 
state  configuration  (V^*0  everywhere)  Is  reached  when  the  flux  tube  be¬ 
comes  sufficiently  kinked  so  that  the  numerator  of  Equation  19  is  zero  for 
each  of  the  three  plasma  elements. 


Assuming  for  the  moment  that  there  are  no  electron  collisions 
(so  that  neo*nei*0),  the  electric  field  Ej_  In  the  quasi-steady  state  limit 
Is  given  by  the  last  term  of  Equation  13,  because  the  first  two  terms 


cancel  In  this  limit.  Figure  4  shows  this  for  the  northern  half  of 

our  model  problem.  Note  that  Ejj^  Is  directed  downward  throughout  the 
central  plasma  element,  but  upward  throughout  the  entire  northern  ele¬ 
ment.  This  allows  the  electrons  to  "ExB"  drift  eastward  In  the  central 
plasma  element,  and  westward  In  the  northern  element,  thus  generating  the 

corresponding  J,  shown  In  Figure  3.  The  Ions  cannot  follow  this  partlcu- 
^  ♦ 
lar  *E*B"  drift  because  the  Ion-neutral  colllslonal  force  v^p^V  In 

Equation  6  opposes  this  E^j^  field  component. 

♦ 

Note  that  the  Ejj^g  vector  shown  In  Figure  4  actually  points  in 
opposite  directions  at  different  positions  along  the  flux  tube,  i.e.,  at 
the  central  and  the  northern  plasma  elements.  Thus  this  quasi-steady 
state  configuration  cannot  be  approximated  by  an  equipotential  solution. 
Furthermore,  38/ 3t  must  be  negligible  in  our  quasi-steady  state,  so  Fara- 
day's  law  requires  VxE*0.  Because  E^j^  is  the  only  electric  field  compo¬ 
nent  transverse-to-B,  this  implies  that  an  E(  must  exist  such  that: 

0  *  7x(Ei+EjJx8)  -  Vx(E  ,+OxB/eN)  (20) 

because  Ejj^  Is  the  last  term  of  Equation  13  (in  the  limit  that  neo  *0). 
This  Is  consistent  with  the  steady-state  assumption  that  V^x  has  no  compo¬ 
nent  In  direction  of  Jx,  I.e.,  that  the  two  terms  on  the  second  line  of 
Equation  12  cancel  one  another  (of  course,  the  first  two  terms  on  the 
r.h.s.  of  Equation  12  also  cancel  one  another  In  the  quasi-steady  state, 
but  not  during  the  initial  transient  stage  when  the  flux  tube  and  plasma 
elements  have  vertical  motions). 

♦ 

The  resultant  quasi -steady  state  E.  fields  are  Indicated  In 
♦  * 

Figure  4.  Note  that  E(  occurs  primarily  at  the  left  and  right  edges  of 
the  plasma  element  where  E,,^.  undergoes  strong  changes.  Furthermore, 
this  E,  exists  across  the  whole  east-west  extent  of  the  plasma  element 


Figure  4.  The  field  and  associated  £„  In  final  quasi-steady  state 
configuration. 
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edges,  and  is  in  opposite  directions  at  the  top  and  bottom  parts  of  the 

elements.  The  E,  is  also  in  opposite  directions  on  the  left  and  right 

sides  of  the  plasma  elements.  Note  that  this  E,  is  not  at  all  correlated 

♦  ■ 

with  the  location  of  J(  ,  and  in  fact  must  exist,  in  this  quasi -steady 
state  configuration,  even  when  electron  collisions  are  negligible. 

Let  us  now  consider  what  would  happen  when  electron  collisions 
are  allowed.  To  keep  close  to  the  Kingfish  data,  suppose  electron  colli¬ 
sions  are  still  negligible  in  the  central  plasma  element  (i.e.,  the  fire¬ 
ball  plasma),  but  extensive  electron  collisions  occur  in  the  region  of  the 
northern  plasma  element  (i.e.,  the  NCR  beta  tube  region).  The  second  term 
of  Equation  12  then  indicates  that  the  flux  tube  velocity  V,-  would  now 

+  -M 

have  an  additional  term  (ng0+ne^ )cJxb/eN.  For  the  northern  plasma  ele¬ 
ment,  this  additional  velocity  is  upward,  and  therefore  allows  the  dis¬ 
torted  field  line  to  relax  toward  the  ambient  (horizontal)  configuration. 
This  motion  seems  consistent  with  the  northward  migration  of  the  Kingfish 
beta  tube  footpoint. 

Electron  collision  effects  in  the  next-to-last  term  of  Equation 

♦  ♦ 

12  indicate  a  further  additional  term  for  equal  to  2(n  /n .  )cJ  /eN. 
Because  ^o*1000neo  »  this  latter  term  is  quite  negligible  in  magnitude 
compared  to  the  previous  additional  term.  Furthermore,  this  term  is 
directed  eastward,  whereas  the  Kingfish  (and  Checkmate)  beta  tube  foot- 
points  clearly  migrated  about  50  km  westward  during  its  northward  motion. 
Thus  electron  collisions  Involved  in  this  term  of  Equation  12  cannot 
explain  the  observed  westward  beta  tube  migration. 

The  observed  westward  migration  would  be  explained  if  the 

special  E  discussed  above  were  to  diminish  in  part  or  wholly.  In  that 
♦  * 

case,  V  would  have  an  additional  term  equal  to  some  fraction  of  -cJ^/eN, 
and  this  Is  correctly  directed  In  the  westward  direction  in  the  northern 
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and  southern  plasma  elements.  The  Kingfish  data  seems  to  imply,  then, 

*► 

that  the  required  special  Ej  may  have  developed  during  the  first  30  sec¬ 
onds  after  the  burst  when  the  beta-tube  stayed  locked  in  place,  but  at 
subsequent  times  it  did  not  reach  the  required  magnitude  so  as  to  cancel 
the  -cJ^/eN  term  in  Equation  12.  The  reason  for  this  is  not  clear,  but 
may  be  related  to  the  substantially  lower  electron  density  N  at  the  beta 
tube  footprints  at  the  later  times. 


SECTION  4 

MAGNETOS PHER I C  PLASMA 


Let  us  now  consider  a  plasma  In  the  magnetosphere  where  colli¬ 
sions  with  neutrals  are  negligible.  We  will  also  neglect  plasma  pressure 
gradients  and  electron  inertia,  so  ion  and  electron  fluid  momentum  balance 
In  the  magnetosphere  Is  approximated  by  the  following  two  simplified 
versions  of  Equations  12-125  and  12-126  of  the  addendum: 


PiVi  »  +  eN(E+V<xB/c)  -  ve^e(V}-Ve) 

(21) 

0  -  -  eN(E+VexB/c)  +  ve1Pe(VrVe) 

(22) 

Addition  of  these  two  equations  yields 

♦ 

•  ♦  ♦ 

PiV1  -  JxB 

(23) 

+ 

This  equation  may  be  solved  for 

*  ♦ 

♦  Pi  ♦  • 

Jx  *  r  bxV1  (24) 


♦ 

In  the  magnetospherlc  case.  It  Is  not  possible  to  eliminate  E.  In  the 

equations  for  Via  and  via  algebra  (In  analogy  with  Equations  10  and  11 

for  the  Ionospheric  case).  Here  we  must  time- Integrate  Equation  23  to 
♦  ♦  ♦ 
obtain  Vu(t),  and  use  Equations  9  and  23  to  relate  V^t)  to  »u(t): 


fU(t) 


V0) 


♦  ♦  , 


/  dt  JxB/p 

A  $ 


(25) 


25 


(26) 


+  V^xb/Q^ 


cJ^/eN 


f-L 


V1J-  +  nei 


cJxb 

eN 


-  -  bx7, 


/(xs,y$,0) 


cJ 


t'uJ_  c  ♦ 

— :  +  i  bxv./E.cbi, 
eN  B  -^1* 


♦  ♦ 
VjxB 


+  n 


H, 

e1  eN 


■¥  ♦ 
JxB 

eN 


(27) 

(28) 


The  last  two  equations  were  obtained  by  solving  Equations  21,  23  and  9 
♦  ♦ 
for  Ej_,  and  using  this  expression  in  Equation  1  to  obtain  Vf4_.  The  above 

Equations  25-28  for  a  magnetospheric  plasma  are  analogous  with  Equations 

10-13  for  an  ionospheric  plasma.  The  two  sets  of  equations  are  quite 

similar,  so  much  of  the  above  discussion  of  flux  tube  motion  adhering  to 

ionospheric  plasma  motion  and  appearance  of  special  E(  fields  also  applies 

to  the  magnetospheric  case. 

For  the  magnetospheric  case  it  is  not  immediately  obvious  that 

“the  electrons  carry  the  current"  (when  one  is  in  an  earth-fixed  frame  of 

reference).  This  does  turn  out  to  be  approximately  the  case  for  simple 

shear  Alfven  waves.  For  example,  Scholer  considered  the  motion  of  a 

6  ♦ 

barium  cloud  In  the  magnetosphere  ,  wherein  he  assumed  that  Jx  was  always 

pointed  In  one  direction  (although  It  could  have  plus  or  minus  signs). 

♦  ♦  ♦ 

Because  Scholer  assumed  V^OJ-O  and  JxB  was  always  in  a  single  direc¬ 
tion,  Equation  25  shows  that  his  did  not  ever  have  a  component  along 
Jx.  Consequently,  one  might  Infer  that  the  electrons  must  carry  the 
current  in  Scholer’s  problem  (in  his  earth-fixed  frame  of  reference).  To 


26 
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actually  prove  this,  one  has  to  retain  the  so-called  “Hall  inertial  term" 
in  the  differential  equations.  The  solution  is  then  very  difficult  to 
obtain,  but  R.  W.  Stagat  has  been  able  to  partly  solve  this  complete 
differential  equation  for  Scholer's  problem. 

Sperling7  has  derived  an  Interesting  analytical  solution  to  an 

extension  of  Scholer's  problem  that  also  includes  a  special  combination  of 

ion  and  electron  collisions  with  background  neutrals.  His  solution 

explicitly  shows  the  deviation  of  the  flux  tube  motion  from  the  plasma 

filament  motion.  In  his  paper,  as  In  Scholer's,  it  was  again  assumed  that 
♦ 

J.  was  always  In  one  direction  (with  plus  or  minus  signs),  so  one  might 

4- 

again  infer  from  Equation  25  that  could  not  have  a  component  along  J^. 

If  one  grants  that  V14.  has  essentially  no  component  along  Jx  in 
the  magnetospherl c  case,  then  Equation  27  shows  that  if  the  flux  tube  is 
to  remain  tied  to  the  plasma  filament,  the  last  two  terms  of  Equation  27 

4- 

must  cancel  each  other.  Consequently,  special  E(  fields  would  also  be 
required  in  the  magnetospheric  plasma,  just  as  in  the  Ionospheric  plasma. 
Also  note  that  the  last  term  of  Equation  28  Is  the  same  E^g  component  of 
Ex  as  discussed  In  the  Ionospheric  case.  Thus  all  the  above  discussion 
for  the  Ionospheric  plasma  can  be  applied  analogously  to  the 
magnetospheric  plasma. 

The  time  evolution  of  the  magnetospheric  plasma  and  magnetic 

♦ 

field  line  configuration  can  be  computed  via  Equations  25  and  27  for  V^x 
and  VfX,  plus  Ampere's  law  and  an  appropriate  equation  for  E(  .  Here 
again.  Equations  26  and  28  are  merely  auxiliary  equations  to  compute  VA. 
and  Ej_  1 f  these  quantities  are  desired. 
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SECTION  5 

SUMMARY  AND  CONCLUSIONS 


For  Ionospheric  plasmas.  It  Is  clear  that  the  electrons  carry 
essentially  all  of  the  transverse-to-B  currents  as  well  as  parallel-to-B 
currents  relative  to  the  neutral  fluid.  For  magnet ospherlc  plasmas  the 
electrons  also  appear  to  be  the  primary  current  carrier,  but  this  case  Is 
not  so  easily  proven  as  that  for  Ionospheric  plasmas. 

If  the  magnetic  flux  tubes  were  to  remain  aligned  with  the 
*► 

moving  Ion  fluid,  special  E(  fields  must  appear  to  cancel  the  curl  of 
EjjxB,  where  this  component  of  Ej_  Is  defined  to  be  the  last  term  of  Equa¬ 
tion  13  or  28: 

*jjxB  s  (29) 


♦  ♦ 

This  special  E(  Is  not  associated  with  J|t  and  It  would  have  to  exist  even 
In  the  absence  of  electron  collisions  If  the  magnetic  flux  tube  and  plasma 
filament  were  to  remain  aligned.  A  consequence  of  these  special  E,  fields 
Is  that  theorists'  predictions  of  Image  strlatlons  In  the  lower  Ionosphere 
would  be  spurious.  Such  Image  strlatlons  have  never  been  observed  In  the 
E-layer. 


Unusual  E,  fields  have  been  observed  above  the  auroral  zone,  and 
i  .  ♦ 

these  fields  seem  similar  to  the  E,  discussed  here. 

The  westward  motion  of  the  Klngflsh  and  Checkmate  beta-tubes, 
during  the  first  three  minutes  after  these  nuclear^  bursts,  may  be  consis¬ 
tent  with  the  Initial  appearance  of  these  special  E,  fields  during  the 
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first  30  seconds  and  their  subsequent  gradual  decay.  The  specific  compu- 

tation  of  the  decay  of  such  E.  fields  is  not  understood,  however,  but  it 

may  involve  the  detailed  evolution  of  the  electron  velocity  distribution 

along  the  whole  flux  tube  in  the  presence  of  the  special  Ej  fields.  On 

the  other  hand,  the  observed  lack  of  initial  westward  motion  during  the 

first  30  seconds  might  also  be  explained  as  due  to  the  relatively  high 

electron  density  at  these  initial  times,  thus  reducing  the  westward 
♦  ♦ 

-cJx/eN  component  of  V^.  in  Equation  12.  At  subsequent  times,  beta-depo¬ 
sition  decreases  and  so  N  also  decreases,  while  Jx  remains  large  at  85  km 

altitude,  thus  allowing  the  beta-tube  foot-points  to  drift  westward.  From 

♦ 

this  point  of  view,  special  E(  fields  would  not  be  required.  A  very 
detailed  3-0  computation  of  the  Kingflsh  and  Checlmate  plasma  evolution 
during  the  first  five  minutes  after  the  bursts  seems  necessary  to  resolve 
whether  or  not  special  E(  fields  are  required  to  explain  the  nuclear  burst 
observations. 

If  the  above  special  Ej  field s  cannot  develop  in  nuclear  burst 

strlatlons,  then  It  seems  likely  that  long  striations  will  gradually  shear 

♦ 

across  adjacent  flux  tubes  in  those  regions  of  the  striation  where  J.  is 
non-zero,  with  a  transverse-to-8  shearing  velocity  equal  to  cJx/eN.  The 
effect  on  the  striated  nuclear  burst  plasma  would  be  to  decrease  the  mean 
electron  density  1?  and  rms  electron  density  fluctuations  ou  at  very  high 
altitudes,  to  increase  the  N  and  at  ionospheric  altitudes,  and  to  ac¬ 
celerate  the  overall  plasma  decay  due  to  the  more  rapid  molecular  chemis¬ 
try  at  ionospheric  heights. 
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APPENDIX 

MOST  SIMPLE  VERSION  OF  IONOSPHERIC  PLASMA  EQUATIONS 


Several  coworkers  have  suggested  to  me  that  It  would  be  more 
♦  ♦  ♦  ♦ 
obvious  and  convincing  that  Is  negligible  compared  to  (l.e., 

that  the  electrons  primarily  carry  the  current  transverse- to-B  relative  to 

♦  ♦  ♦ 

the  neutrals)  if  this  were  proved  using  the  standard  E.  and  E.xb  basis 
vectors  for  V.  ,  V  ,  and  J  .  Let  us  therefore  use  a  local  reference  frame 
moving  with  the  neutral  fluid  so  Vo*0.  We  also  neglect  electron  colli¬ 
sions  (neo*0  and  ne<-0).  Ion  Inertia,  ^(P^+Pg),  g,  etc.  Keeping  only 
Ion-neutral  collisions  (n^0),  Equations  12-133,  12-134,  and  12-135  yield 
the  standard  formulas  for  this  most  simple  Ionospheric  plasma  motion: 
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Taking  the  dot  product  of  Jx  with  and  ,  we  then  find: 
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Therefore  only  the  electrons  move  In  direction  of  Jj_  relative  to 
the  neutrals.  The  ions  do  not  move  In  direction  of  J.  ,  but  do  move  In 
direction  of  Jjxb  relative  to  the  neutrals.  This  is  consistent  with 
Equations  10  and  11  of  Section  2;  those  equations  are  more  general  in  that 
they  also  Include  the  effect  of  electron  collisions,  but  since 
Hi  "100<Vj  ,  It  Is  clear  that  the  electrons  are  the  primary  current  car¬ 

rier  transverse-to-B  as  well  as  parallel-to-B. 

Figure  1-1  shows  a  pictorial  version  of  Equations  1-1.  1-2.  and 

1-3  suggested  to  the  author  by  C.  Longmire.  We  have  drawn  the  vector 

diagrams  for  various  nio  values  (1.e.»  for  various  Ionospheric  altitudes) 

and  adjusted  the  E^  magnitude  and  direction  so  that  .vts  the  same  at  all 

altitudes  (l.e.,  for  all  n^Q  values).  Note  that  the  Jx  vector  Is  always 

perpendicular  to  and  points  In  the  same  direction  at  all  altitudes  In 

this  simple  example  wherein  Vq«0.  Thus  in  this  simple  case  the  Ion  column 

moves  as  a  unit  transverse-to-B,  but  the  electric  field  cannot  be  taken 

♦ 

as  an  equl potential  because  E^  rotates  and  changes  Its  magnitude  at  the 
various  altitudes  (l.e.,  various  n^Q  values). 


i 

i 


Figure  1-1.  Relationship  between  Vu,  v^,  cE^/B,  and  cJx/eN  at  various 
altitudes.  Vu  has  been  chosen  to  have  the  same  direction  and  magnitude 
at  all  altitudes;  this  then  requires  that  Ex  rotates  clockwise  and 
Increases  In  magnitude  (as  /l+nJQ)  as  one  descends  In  altitude.  Note 
that  Is  always  perpendicular  to  Vix<  Here  It  has  been  assumed  that 
VO  at  all  altitudes. 
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ADDENDUM:  Pages  617-621  from  DNA4501F 


12.4.4  LOW-6  PLASMAS  ABOVE  100  KM  ALTITUDE 

In  the  normal  ionospheric  E  and  F-layers  the  ion-electron  pressure  P  is 
much  less  than  the  magnetic  pressure  B2/8ir.  The  ratio  8ttP/B2  is  often 
referred  to  as  the  plasma-B.  Such  low- 8  plasmas  also  occur  after  a  few 
minutes  subsequent  to  high  altitude  nuclear  bursts  and  for  barium  cloud 
releases  in  the  ionosphere.  Under  these  conditions  it  is  not  possible  ener¬ 
getically  to  cause  substantial  distortions  of  the  magnetic  field.  Thus  3$/3t 
is  very  small,  and. as  a  consequence  of  Maxwell's  Equation  (12-2),  the  elec¬ 
tric  field  must  be  primarily  electrostatic.  Experience  with  numerical  com¬ 
putations  shows  that  it  is  very  difficult  to  accurately  solve  the  magneto¬ 
hydrodynamic  Equations  (12-100)  through  (12-107)  when  the  plasma-8  is 
small.  One  must  recast  these  equations  to  obtain  a  more  viable  procedure 
for  finding  a  solution  (Ref.  12-4) . 

A  key  point  is  that  currents  will  be  induced  in  the  high  altitude  plasma  as 
it  moves  about  in  response  to  the  applied  forces  such  as  gravity,  neutral 
winds,  pressure  gradients,  etc.  The  effect  of  these  currents  is  to  average 
the  applied  forces  along  each  field  line  via  J*B  forces  such  that 
essentially  all  of  the  electrons  in  a  given  magnetic  flux  tube  will  E*B 
drift  together  to  another  field-aligned  flux  tube.  Because  only  a  negligibly 
small  space-charge  is  needed  to  set  up  the  spatially  varying  electrostatic 
field  I  that  causes  the  1=x1j  drift  of  the  plasma,  it  is  crucial  that  the 
current  density  5  be  essentially  divergence-free,  i.e.,  that: 

V-3  *  0  (12-123) 

To  employ  this  equation  we  must  express  3  in  tents  of  the  electric  field 
and  the  forces  acting  on  the  plasma.  For  the  conditions  discussed  in  this 
subsection,  the  ion-electron  plasma  is  collisional  because  the  plasma 
evolves  on  a  time  scale  of  minutes  and  the  expected  plasma  temperature- 
density  parameters  are  in  the  range  of  less  than  0.1  eV  at  104cm~s  to 
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less  than  2  eV  at  109cm~3.  On  the  other  hand,  above  100  km  altitude, 
collisions  with  neutrals  are  insufficient  to  provide  tight  coupling  with 
the  neutral  fluid.  Thus  we  must  employ  a  multi-fluid  description,  with 
the  neutrals,  the  ions,  and  the  electrons  obeying  the  following  momentum 
conservation  equations: 
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(12-125) 


(12-126) 


The  subscripts  o,  i,  and  e  respectively  indicate  neutral,  ion,  and 
electron  fluid  parameters.  The  ion-neutral,  electron-neutral,  and  electron- 
ion  collision  frequencies  are  denoted  by  v.  ,  v  ,  and  v  ..  Detailed 

XO  00  ox 

calculations  of  these  collision  frequencies  for  nuclear  burst  conditions 
are  given  in  Ref.  12-5;  approximate  values  are  indicated  in  the  glossary 
of  Chapter  13  of  this  compendium. 


The  terms  involving  v^0»ve0»  and  ve^  in  the  three  above  equations 
account  for  momentum  transfer  between  the  three  fluids  via  collisions. 

In  Equation  (12-124),  the  last  term  is  negligible  compared  to  the  others. 
The  usual  magnetohydrodynamic  equation  is  obtained  by  adding  Equations 
(12-125)  and  (12-126).  Neglecting  some  small  terms,  this  yields: 
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(12-127) 
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Equations  (12-124)  and  (12-127),  together  with  (12-106)  and  (12-107),  are 
referred  to  as  the  two-fluid  MHD  equations. 

Returning  now  to  the  problem  of  obtaining  a  suitable  expression  for  J,  we 

will  solve  Equations  (12-125)  and  (12-126)  for  1.  and  V  and  then  find 

i  e 

J  from 


j  *  f  (VV 


(12-128) 


To  solve  those  two  equations  for  the  transverse-to-1  velocity  components 

V.  and  1  we  introduce  the  auxiliary  notation: 
ix  ex 
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(12-131) 


d  s  (l^0)(i^o)  ♦  nei(nio^eo)  [>2T1.oneoxnei(nioxneo)l  (12-132) 


where  the  n's  are  the  ratios  of  the  collision  frequencies  to  the  appro¬ 
priate  gyro-frequency.  These  n's  are  all  positive  numbers.  For  nuclear 
burst  conditions,  n  and  n4_  are  less  than  unity  above  about  70  km 
and  130  km  altitude  respectively,  and  they  are  nearly  proportional  to  the 
neutral  fluid  density.  At  any  particular  point  the  ratio  ne0/n^0  is 
small  and  about  equal  to  1/1000.  The  nei  is  proportional  to  N#  and 
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at  altitudes 


outside  fireballs  the  n  *  is  generally  smaller  than  ri 

61  60 

below  200  km.  Above  300  km  altitude,  or  inside  fireballs,  ri  .  is  usually 

61 

larger  than  n  and  may  also  be  larger  than  n.  .  At  altitudes  above 
eo  10 

about  ISO  km  d  is  approximately  unity  because  n-Q»  n  Q>  and  ngi  are 
generally  small  at  high  altitudes  . 

The  exact  solutions  for  V.  and  V  are  obtained  by  inverting  a  4x4 

li  ei  ° 

matrix,  with  the  result  that 
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in  Equation  (12-128),  we  obtain: 


1 


Using  these  expressions  for  and 
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(12-135) 


Equations  (12-125)  and  (12-126)  may  also  be  solved  for  the  parallel-to-B 
current  density  J  : 

®  [%»  *  .*  £ )] J.  ’  f  ft  ♦**.)-£  * 


(12-136) 


No  approximations  have  been  made  in  deriving  Equations  (12-133) -(12-136) 
from  the  MID  Equations  (12-125)  and  (12-126)  except  the  negligible  one  that 
N.»N  »N.  In  particular,  the  £  fields  appearing  in  these  equations  could 
be  inductive  as  well  as  electrostatic. 
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